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1 . INTRODUCTION 


Infrared  radiation  from  clouds  is  an  important  background  signal 
contribution  that  must  be  considered  in  the  design  of  high-altitude  sensor 
systems  deployed  for  missile  and  aircraft  detection.  Thermal  radiation 
from  clouds  and  its  transmittance  through  the  atmosphere  can  be  computed 
in  a relatively  straightforward  manner.  On  the  other  hand,  the  calculation 
of  apparent  cloud  radiance  that  results  from  the  diffuse  reflection  of  solar 
radiation  from  clouds  is  a much  more  complicated  problem  in  multiple  - 
scattering  radiative  transfer  theory.  Even  with  the  simplifying  assump- 
tions of  an  infinite,  plane -parallel  cloud  layer  composed  of  spherical 
particles,  the  computation  of  exact  solutions  of  the  reflected  radiation 
field  for  arbitrary  solar  and  viewing  angles  is  complex  and  time  consum- 
ing, prohibitively  so  for  routine  application  to  parametric  system  design 
studies.  The  overall  objective  of  the  present  study  is  to  investigate 
approximate  cloud  reflectance  models  that  can  be  effectively  applied  in 
systems  analysis. 

The  accuracy  of  these  models  is  not  of  paramount  importance.  The 
problem  is  usually  phrased  by  the  system  analyst  as  "how  much  reflectance 
would  I expect  from  high-altitude  clouds"  without  any  specification  of  cloud 
type,  homogeneity  of  cloud  cover,  particle  type  or  size  distribution,  etc. 

At  the  same  time,  very  large  errors  cannot  be  acceptable  since  these 
could  presumably  influence  the  system  design.  The  subjective  criterion 
set  for  acceptable  model  accuracy  here  is  that  the  model  should  predict 
within  a factor  of  two  of  the  exact  reflectance  for  a typical  cloud  type.  The 
cloud  type  chosen  is  a model,  plane -parallel,  semi-infinite  cumulus  cloud 
composed  of  spherical  water  or  ice  particles.  Further  consideration  of 
this  test  standard  is  made  in  Section  3. 


The  method  of  approach  to  the  problem  is  straightforward.  A liter- 
ature search  was  made  for  suitable  models  that  appeared  amenable  to  the 
problem.  The  methods  were  adapted  to  the  present  problem  in  the  simplest 
manner  possible  and  the  predictions  for  the  test  standard  compared  to  the 
exact  results  for  the  standard  case. 


2ir  1 

L(t,P,0)  ^ ^ p(p,0,p',0')L(T,p«,0«)dp'd0' 

6 -i 


Note  that  the  subscript  D has  been  dropped.  From  this  point  on,  the  sym- 
bol L will  be  used  to  denote  L^.  The  function  p(p,  0,p' ,0')  is  the  scattering 
phase  function  for  the  cloud  particles  and  describes  the  probability  for 
scattering  by  a single  particle  from  the  direction  p'0'  to  the  direction  p0. 
For  spherical  particles,  p does  not  depend  individually  on  p and  0 but  is  a 
function  only  of  the  scattering  angle  cp  defined  by 

I 2 2 ' 

cos<P  = cos(0-0')\ /(1-p  )( 1-p'  )+pp'.  (2) 


- x-Erte 

4ir  0 


p(p,0,  “M'q#  0q)* 


2.  EQUATION  OF  RADIATIVE  TRANSFER 


Consider  parallel  solar  radiation  of  spectral  flux  density  E^W/cm^-pm) 
impinging  on  the  top  of  a plane -parallel,  homogeneous  cloud  layer  at  zenith 
angle  0q  (see  Fig.  1).^*^*^  The  incoming  radiation  is  scattered  by  the 
particles  comprising  the  cloud,  and  a steady  state  spectral  radiance  field 
MT»  p*  0)(W7cm^-sr-pm)  is  established  within  the  cloud,  t is  a dimension- 
less measure  of  distance  into  the  cloud  (from  the  top)  and  is  defined  later, 
p is  the  zenith  angle  variable  p = cos6.  The  field  L is  composed  of  two 
parts,  a direct  component  Lq  which  contributes  to  L only  for  0 = ir  - 0q 
and  0 = 0q,  and  a diffuse  component  L^.  The  solution  for  the  direct  com- 
ponent is  Lq  = EQexp(-pT).  The  diffuse  component  for  a plane -parallel, 
vertically  homogeneous  cloud  of  infinite  horizontal  extent  under  conditions 
of  thermodynamic  equilibrium  and  with  no  embedded  sources  of  radiation 
is  governed  by  the  following  integro-differential  equation  of  radiative 
transfer: 


MIINMMmi 


The  phase  function  is  normalized  such  that 


2ir  1 


[ ^ p(p,0,p',0')dp'd0'  = 2ir  J p(x)dx  = 4nu)Q. 

o -i  -l 


(3) 


where  x = cosep.  uu^  is  the  single-scattering  albedo  of  the  cloud  particles 
and  is  defined  in  terms  of  the  particle  cross  sections  for  absorption  ( a ) 
and  scattering  (a#)  by 


<u. 


(4) 


where  o is  the  total  extinction  cross  section  a -a  + a . These  cross 
e e s a 

sections,  and  the  differential  scattering  cross  section  da#(«p)/dQwhich 
defines  the  phase  function  through 

4_  d°(<P) 

p(cosCp)  = - gjj-  (5) 

e 

are  considered  in  the  following  section.  The  optical  depth  parameter  t is 
defined  by  t = Noez  where  N is  the  cloud  particle  density. 

The  boundary  conditions  that  the  radiance  field  L must  satisfy  are 
(1)  that  there  be  no  downward  directed  diffuse  radiance  at  the  top  of  the 
cloud,  i.  e.  , 


L(0,  p,  0)  - 0 p<0,  (6! 

and  (2)  for  a finite  thickness  cloud,  that  there  be  no  upward  directed 
radiance  at  the  bottom,  i.  e.  , 

L(t0.h,*)  = 0 ji>0,  (7; 


or  (3),  for  a semi- infinite  cloud,  that  L(t)  remain  bounded  as  t - 

The  three  terms  on  the  right  of  Eq.  (1)  have  clear  physical  signifi- 
cance. The  first  term  reflects  Beer's  law  that  the  change  in  L is  propor- 
tional to  L.  The  second  (integral)  term  is  a source  term  corresponding 
to  the  scattering  of  diffuse  radiation  from  all  directions  into  the  particular 


— »_« 


direction  p0,  and  the  last  term  is  a source  function  corresponding  to  the 
scattering  of  radiation  from  the  direct  radiance  component  Lq  into  the 
direction  p0. 

Given  the  phase  function  p(cos<p),  its  normalization  constant  u)q  and 
the  boundary  conditions  of  Eqs.  (6)  and  (7),  the  scattering  problem 
described  by  Eq.  (1)  is  fully  defined.  In  principle,  one  solves  this  equa- 
tion for  L(T,p,0).  Then,  an  appropriate  cloud  reflection  coefficient 
p(sr  * ) is  defined  by 


p > 0. 


In  practice,  the  solution  of  Eq.  (1)  is  formidable.  The  primary  purpose 
of  this  paper  is  to  investigate  rational  approximate  solutions.  Descriptions 
of  the  approximations  comprise  Section  4 of  this  report.  For  the  most 
part,  the  difficulty  in  solving  Eq.  (1)  is  caused  by  the  presence  of  the  inte- 
gral scattering  source  term,  and  the  principle  methods  of  approximation 
simplify  this  term.  The  first-order  approximation  of  completely  neglect- 
ing this  term  is  the  so-called  single -scattering  approximation  and  is  also 
considered. 


3.  TEST  CONDITIONS 


The  cloud  standard  is  a plane -parallel,  semi-infinite  cumulus  cloud 
composed  of  spherical  water  droplets.  The  size  distribution  of  particle 
radii  is  described  by^'  ^ 

y 

n(r)  = ar°e"^r  (9) 

where 

a = N{vP("+1,/v/r(2ii)]. 


N = particle  density  (cm-3),  rc  = radius  of  peak  distribution  (pm),  and  a 
and  y are  size  parameters  controlling  the  shape  (e.  g.  , width  and  skew- 
ness) of  the  distribution.  For  the  present  work,  a = 6,  y = 1 and  r^.  = 4pm 
are  used  and  the  resulting  distribution  f(r)  = n(r)/N  is  shown  in  Fig.  2. 

Necessary  inputs  to  any  nonempirical  model  of  cloud  reflectance  are 
the  scattering  (total  and  differential)  and  absorption  cross  sections  for 
single  cloud  particles.  These  parameters  were  computed  according  to 
standard  Mie  scattering  theory'  ’ and  averaged  over  the  size  distribution 
of  Eq.  (1).  The  index  of  refraction  for  water  was  taken  from  Ref.  (6)  and 
is  shown  in  Fig.  3.  Results  in  the  2-20pm  spectral  region  for  the  total 
scattering  (a  ),  total  absorption  (a  ) and  total  extinction  (a  - a + a ) 
cross  sections  are  shown  in  Fig.  4.  The  differential  scattering  cross 
section  dcs(<p)/dO  as  a function  of  scattering  angle  at  \ = 2.  70  and  4.  35pm 
is  shown  in  Fig.  5.  A fundamental  feature  of  the  differential  scattering 
cross  sections  (i.  e.,  the  phase  function)  evident  in  Fig.  5 is  the  strong 
probability  for  forward  scattering.  Advantage  is  taken  of  this  feature  in 

*S.  J.  Young,  Calculation  of  Scattering  Cross  Sections  for  Clouds,  ATM 
75(5409-40)-5,  The  Aerospace  Corporation,  El  Segundo,  California,  9 
June  1975. 
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the  formulation  of  various  strong-forward-scattering  approximations  in  the 
next  section.  The  principle  index  for  measuring  the  degree  of  forward 
scattering  is  the  fraction  of  scattering  events  for  which  the  scattering  angle 
is  less  than  90  degrees.  The  fraction  T)  is  shown  in  Fig.  6.  Also,  the  type 
of  simplifications  that  can  be  made  in  modeling  cloud  reflectance  is  strongly 
influenced  by  the  albedo  for  single  particle  scattering,  a parameter  that 
measures  the  fraction  of  events  in  which  the  incident  radiation  is  scattered 
rather  than  absorbed.  This  parameter  is  defined  by  <i>q  = o^lo^  and  is 
shown  in  Fig.  7. 

Exact  calculations  for  reflectance  from  this  model  cloud  have  been 
made  by  Hansen  etal.^"*^  using  the  "doubling"  method  of  radiative  trans- 
fer. The  test  case  against  which  the  results  of  the  approximate  methods 
are  compared  is  shown  in  Fig.  8.  The  case  is  for  normally  incident  radia- 
tion at  X = 3.  4 um. 


4.  MODEL  DESCRIPTIONS  AND  RESULTS 


I 


A.  Single -Scattering  Model 

The  single-scattering  model^^  for  cloud  reflectance  considers  the 
diffuse  radiance  field  within  the  cloud  to  be  composed  only  of  those  photons 
that  have  been  scattered  once  from  the  primary  incident  beam.  The  radia- 
tive transfer  equation  for  this  condition  is  obtained  from  Eq.  (1)  by 
neglecting' the  integral  source  scattering  term.  The  solution  for  L is 
obtained  simply  by  solving  the  resulting  first-order  differential  equation. 
The  result  for  a semi-infinite  cloud  that  satisfies  the  boundary  conditions 
MO,  p,  0)  = 0 for  p < 0 and  L(T,p.,  0)  -»  0 for  t -♦  ® is 


L(t,h,0)  = 


TT  P(cos<P0)[e-T/>i0  . eT/**] 


'0  


^p(coscp0)e 


(JL  < 0 


> 0 


(10) 


where 


cos<p0  = cos(0-0Q) 


)(l-p0)  - fipQ. 


The  reflection  coefficient  of  Eq.  (8)  is  then 


(ID 


p(^0)  = 4^p(*os  Virnv  (12) 

The  comparison  of  the  prediction  of  this  result  with  the  standard 
test  case  is  shown  in  Fig.  9.  The  single -scattering  approximation  under- 
predicts the  reflectance  by  about  a factor  of  three. 

Although  the  single -scattering  approximation  fails  as  a general 
approximation,  it  has  two  redeeming  features;  (1)  it  is  extremely  simple 
to  apply,  and  (2)  it  can  be  quite  accurate  in  some  spectral  regions.  The 
first  feature  is  used  to  advantage  in  the  formulation  of  modified  single- 
scattering approximations  (Sections  B and  C).  Results  from  a previous 


SINGLE  SCATTERING 


Fig.  9.  Comparison  of  Single -Scattering  Approximations 
with  Test  Case 


study  on  cloud  reflectance  has  established  that  this  model  is  adequate  for 
practical  application  in  the  2-  to  3- pm  spectral  region.  In  Fig.  10,  a 
comparison  is  made  between  the  model  and  experimental  field  measure- 
ments. ^ ^ The  measurements  were  made  with  an  airborne  sensor 

flying  1 km  above  a uniform  cloud  deck  with  an  observation  zenith  angle 
of  80  degrees,  a solar  zenith  angle  of  88.  5 degrees,  and  a nearly  coplanar 
geometry  for  the  entrance  and  exit  rays  (i.  e. , 0 - 0q  “*  0).  When  proper 
account  of  atmospheric  attenuation  is  made  of  the  path  from  space,  to  the 
cloud,  and  back  to  the  sensor,  the  agreement  is  quite  good. 

The  success  of  the  single -scattering  model  in  the  2.  7-pm  region  can 
be  attributed  to  the  large  value  of  the  absorption  coefficient  of  liquid  water 
in  this  spectral  region  (Fig.  3b).  This  large  value  is  carried  over  into  a 
broad  peak  in  the  spectrum  for  the  absorption  cross  section  oa  (Fig.  4) 
and,  consequently,  into  a broad  valley  in  the  spectrum  for  the  single  - 
scattering  albedo  u>q  (Fig.  7).  The  effect  is  that  photons  that  must  be 
scattered  more  than  once  in  order  to  escape  the  top  of  the  cloud  in  the 
desired  reflection  direction  have  a large  probability  of  being  absorbed 
rather  than  scattered.  At  \.  = 3.  4-pm,  on  the  other  hand,  the  absorptance 
by  liquid  water  is  small,  and  higher  orders  of  scattering  contribute  sig- 
nificantly to  the  reflected  radiation  field. 

B.  Two-Stream  Theoretical  Modification  of  the  Single -Scattering  Model 


The  general  development  of  radiative  transfer  in  scattering  media 
has  produced  a set  of  approximate  solutions  collectively  known  as  "flux 
approximations.  " Among  the  more  important  are  those  due  to  Eddington, 
Schuster  and  Schwarzchild  (two-stream  approximation),  and  Sagan  and 
Pollack  (modified  two-stream  approximation).  Use  of  these  approxima- 
tions for  cloud  reflectance  has  been  considered  by  Irvine.  The  essence 

of  these  approximations  is  the  sacrifice  of  angular  information  of  the 

*S.  J.  Young,  Solar  Scattering  Calculations,  ATM-75(5409-41)-3,  The 
Aerospace  Corporation,  El  Segundo,  California,  14  December  1974. 
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Fig.  10.  Comparison  between  Single -Scattering  Model  and 
Experimental  Data  in  the  2.  7 -pm  Spectral  Region 


general  radiation  field  in  order  to  efficiently  compute  the  total  flux  in  a 
given  direction.  In  application  to  the  cloud  reflectance  problem,  these 
approximations  yield  only  the  albedo  of  the  cloud,  that  is,  only  the  total 
radiation  reflected  into  the  upper  hemisphere.  Because  of  this  lack  of 
angular  information,  these  types  of  approximations  are  not  considered 
here  as  final  models.  This  approach  is  used,  however,  in  order  to  derive 
an  improvement  on  the  simple  single -scattering  model. 

The  basic  simplification  made  in  Eq.  (1)  to  obtain  these  flux  approx- 
imations is  to  assume  a phase  function  of  the  form 

p(ti,  0,ti',0')  = 41TUUQ1  r)6  (|xf— m.)6(0'-0)  + (l-i’|)6(p.+ti>)6(0l-ir-0)  • (13) 

where  6 is  the  Dirac  delta  function.  This  form  allows  only  direct  forward 
and  backscatter  along  the  primed  beam  axis,  ti  is  the  forward  scattering 
parameter.  The  radiance  field  is  similarly  constructed  as 

L(t,h,0)  = rE+(T)6(p-po)6(0-ir-0o)  + E’(t)6(p+po)6(0-0o)]  . (14) 


Thus,  the  field  is  assumed  to  contain  components  only  in  and  opposed  to 
the  direction  of  the  incident  solar  radiation.  E+  and  E are  the  radiance 
flux  variables  for  the  forward  and  backward  directions,  respectively. 
Substitution  of  Eq.  (13)  and  (14)  into  Eq.  (1)  yields  the  following  pair  of 
coupled,  first-order  differential  equations  governing  one-dimensional 
radiative  transfer 


dE^li  = 

dT 

. <*E>) 


E"(t) 


u>0(  1 -rj)' 

. *0  . 

U>0(l-T|) 

■ r^r 


E+(t)  - u,0E0(l-ri)e_T/tA0 
■]e+(t)  - u,0E0rie"T/»10  . 


(15) 


The  relevant  boundary  conditions  are 

E~(0)  = 0 

and 

E*(t)  - 0 as  t - ® . 
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The  solution  of  this  set  of  equations  yields* 

EV)  = „ r e-kT/^o 

Eq  *Voe 

£'lTi  = „ re-kT/^0  . e"T/M’o"l 

EQ  ^Oi  C 6 j 


whe  re 


k = / (l-ou0)C  1-(2t|-1  )0JqD  , 


k - (1-Wq) 
T0  = k + (1-u)q)  * 


The  reflection  coefficient  (total  albedo)  is  obtained  as 


E (0) 

P = — ^ — L = u.r.. 
^ Eq  0 0 


Adamson  1 ' has  expanded  Tq  in  orders  of  scattering  to  obtain 


'o  = E Pn 


where 


[n/2] 


■ -8®“^  E [-  • 

j=0  v 

The  contributions  to  the  albedo  from  once,  twice,  and  three  times  scattered 
photons  are,  respectively 
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Pj  * 2 

*i  u>2 

p2  =^^2n2(i-Ti), 
2 


P3  = — 5-  [ (2ri)  2(1  -n)  + 4 (1-ri)  . 

2 

In  the  limit  of  strong  forward  scattering  (r|  -•  1),  these  contributions 
approach 

•>1— 

p2— 


and,  in  general 


•HoW 


D_u,n  (i.  ml 

Pn  P’0  0 2 


In  this  limit,  one  can  calculate  simply  the  ratio  of  all  multiple  scattering 
contributions  to  the  first  order  (single-scattering)  contribution  as 


C = — 


X pn 


lHB0  * 


Although  this  result  has  been  obtained  in  the  flux  approximation,  it  is 
applied  as  a correction  factor  to  the  single -scattering  result  of  Eq.  (12) 
to  obtain  an  improved  model 
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p(H»0)  = Cpgg(ji,0)  (27) 

where  pgg  is  the  single-scattering  approximation.  The  form  of  the  cor- 
rection factor  is  intuitively  correct.  As  the  probability  for  single  scattering 
increases,  i.  e. , as  uu^  increases  towards  unity,  the  effects  of  multiple 
scattering  contributions  to  the  reflectance  must  increase.  This  intuitive 
result  is  evident  in  the  increase  of  the  multiple-scattering  correction 
factor  C with  u)q.  This  approximation  must  fail,  however,  for  conserva- 
tive scattering  (ui^  -•  1)  since  then  p -• «.  For  the  present  application,  this 
poses  no  problem  since  u)q  never  exceeds  ~ 0.  9 (see  Fig.  7). 

For  the  test  case  condition  of  Fig.  8,  u>q  = 0.729  and  the  correction 
factor  is  C = 3.  56.  The  product  of  C and  pgg  is  shown  in  Fig.  9.  The 
agreement  with  the  exact  reflectance  curve  has  been  improved  so  that  the 
maximum  error,  except  at  0 *■  90  degrees,  is  only  ~ 50  percent. 

C.  Empirical  Modification  of  the  Single -Scattering  Model 

Comparison  of  the  predictions  of  the  theoretically  modified  single - 
scattering  model  with  exact  calculations  was  extended  beyond  the  single 
comparison  of  Fig.  9 in  order  to  determine  the  range  of  its  application. 

The  exact  results  were  obtained  from  Refs.  9 and  11.  The  results  of  this 
comparison  are  given  in  Table  1 and  Fig.  11.  Column  one  of  Table  1 gives 
the  wavelength  of  the  eight  cases  used  in  the  comparison.  The  second 
column  lists  the  figure  number  and  reference  from  which  the  exact  result 
was  obtained.  The  size  distribution  function  of  Eq.  (9)  with  r^.  = 4 pm, 

V = 1 and  or  = 6 applies  in  all  cases,  but  for  three  of  the  cases,  the  cloud 
particles  were  assumed  to  be  spherical  ice  particles.  The  cloud  particle 
type  is  listed  in  column  three  of  the  table.  The  fourth  column  tabulates 
the  single-scattering  albedo  for  the  wavelength  X.,  the  fifth  column  is  the 
theoretical  correction  factor  of  Eq.  (26),  and  the  last  column  is  the  factor 
obtained  from  the  exact  results.  This  last  factor  was  obtained  from  the 
figures  referenced  in  the  second  column  by  computing  the  mean  ratio 
between  the  exact  and  single-scattering  reflectance  over  whatever  variable 
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Table  1.  Data  for  Construction  of  Empirically  Modified 
Single  *Scatte ring  Model. 


\(pm)  Fig.  /Ref. 


Cloud 

Type 


Exact 


d/<i-.0i  ins,; 


r 


I 


I 


(e.  g. , 9 or  0)  was  used  to  plot  the  results  in  Refs.  9 and  11.  A plot  of 
this  exact  factor  versus  the  theoretical  factor  1 /(I -u)Q)  is  shown  in  Fig.  11. 
For  u)Q  < 0.  9[  1/(1  -iu0)  < 10],  the  comparison  between  the  model  and  exact 
results  is  very  good.  For  more  conservative  scattering,  however,  the 
model  overpredicts  the  multiple -scattering  contribution  to  reflectance. 

The  model  predicts  a continuing  linear  increase  with  1/(1  -«Jq)  while  the 
exact  results  indicate  a saturation  for  > 0.  99[  1 /(I -U)q)  ^ 100]  such 
that  the  true  reflectance  is  never  more  than  25  times  larger  than  the 
single -scattering  prediction.  An  empirically  modified  single -scattering 
model  was  thus  constructed  by  taking  the  correction  factor  to  be  its 
theoretical  value  for  uuQ  £ 0.  9,  the  saturation  value  Cg  = 25  for  S 0.  984, 
and  the  interpolated  value 

C = 10F  (28) 

e 

with 

F = 1.  3979  - 0.6282  riog10(^^)]  (28) 

for  0.  900  £ ou-  £ 0.  984.  This  correlation  for  C is  the  curve  labeled 
0 ® 

Empirical  Factor  in  Fig.  11.  In  all,  the  approximation  appears  to  be 
accurate  to  ^ ± 50  percent  for  most  cases,  with  an  occasional  excursion 
to  a factor  of  two  error. 

An  example  of  the  application  in  the  2-  to  20-pm  region  of  the  single- 
scattering, theoretically  modified  single -scattering  and  empirically  modi- 
fied s ingle -scattering  models  is  shown  in  Fig.  12.  The  sun  is  at  the  zenith 
and  the  observation  is  at  the  coplanar  45  degrees  zenith. 

D.  Turner  Model 

Turner^’  has  developed  and  applied  a multiple -scattering 

model  that  is  developed  along  the  lines  of  the  flux  approximations,  but 
which  retains  the  dependence  of  the  reflectance  on  the  entrance  and  exit 
zenith  angles.  In  effect,  this  model  is  the  next  logical  step  after  the  flux 


A 
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approximations.  The  solutions  for  E±(t)  of  Eq.  (17)  are  substituted  into 
Eq.  (14)  for  L(t,|i,  0)  and  this  result,  in  turn,  substituted  into  the  radiative 
transfer  equation  [Eq.  (1)].  Equation  (13)  for  p(}i,  0,  p1 , 0' ) is  also  sub- 
stituted into  Eq.  (1),  but  this  time,  only  into  the  integral  scattering  term. 
Integration  over  the  angular  variables  in  the  integral  scattering  term 
reduces  Eq.  (1)  to  a linear  first-order  differential  equation  for  which  the 
solution  is  (for  a semi-infinite  cloud) 


4^rP(cosV  + r0P(-co«v]^^Tk?/c"kT/>i0  " eTAl]  u<0 

(2 

4F  i P<coflV  + roP(-Coa<pO)]p0^kpe'kT/^°  u > 0 


from  which  the  reflection  coefficient  is  obtained  as 


1 T *1  Mn 

p(P*0)  = 4^  P^V  + r0P<-co8V.  nQ  + kp  ti>  0 (30) 

where  cpQ  is  defined  by  Eq.  (11),  k by  Eq.  (18),  and  rQ  by  Eq.  (19). 

The  significant  feature  of  this  model  is  that  it  augments  the  contri- 
bution from  single -scattering  by  adding  a "backscatter  component,"  that 
is,  the  component  r^pt-coscpp).  The  prediction  of  this  model  for  the 
standard  test  case  is  shown  in  Fig.  13.  For  this  normal  incident  case, 
the  principle  failing  of  the  model  is  evident  — a gross  overprediction  of 
the  reflectance  results  for  direct  back  reflection  (0  = 0).  This  results 
because  the  Turner  model's  "backscatter"  contribution  to  the  single- 
scattering contribution  is,  in  fact,  the  strong  forward  scattering  lobe  of 
the  phase  function.  For  large  exit  zenith  angles,  this  model  adds  little 
to  the  predictions  of  the  single-scattering  model  (compare  Figs.  9 and  13 
for  0 *■  90  degrees). 
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E.  Romanova  Approximation 

For  scattering  problems  involving  phase  functions  that  are  strongly 
peaked  in  the  forward  direction,  Romanova^ ^ 2*^  gives  a procedure  for 
separating  the  general  transfer  equation  into  two  parts.  One  part  describes 
the  strongly  angle-dependent  forward  radiation  field,  while  the  other 
describes  the  weakly  angle -dependent  remainder  field.  The  separation  is 
effected  by  rewriting  the  transfer  equation  [Eq.  (1)]  as 


2ir  1 

-ti0  ^ p(H,0,p',0')L(T,p',0')dp'd0' 

b -i 


(31) 


- T?e-T/‘i°P(^^-^o*0o)  ■ c(ixo  + d1^"1  • 


When  the  parameter  e is  set  to  unity,  the  exact  transfer  equation  is 
recovered.  When  e = 0,  the  equation  describes  the  field  in  the  direction 
p.  “•  — Pq.  The  field  variable  L(T,p,  0)  is  now  decomposed  by 

L(t,  h,  0)  = Lj  (r,  fi,  0)  + «L2(t,  ft , 0)  (32) 


and  substituted  into  Eq.  (31).  Collecting  terms  independent  of  e yields 
the  following  transfer  equation  for  Lj 

dL  (T,fi,0)  . 2p 

^ = Lj(T,p,0)  - ^ \ \ p(fl,0,fl',0')L(T,fl',0')dHL'd0' 

0 -1 


■^0 


E0  _-t/k-0 


(33) 


4ir 


p(fi,  0,  -HO,0O)- 


Collecting  terms  containing  e and  setting  e = 1 yields  the  following  transfer 
equation  for  L ^ 

dL2(T,fi,0)  2p  i 

H = L,(T,fji,0)  - ^ \ \ p(p, 0,  p' , 0' )L?(t, p' , 0' )d(i'd0' 


0 -1 


dL.(T,|JL,0) 

-(^o  + I*) 37 


(34) 
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Presumably,  Eq.  (33)  describes  the  field  in  the  direction  p •*  -p^  and 
0 *•  0q  (which  is  strongly  angle  dependent,  at  least  for  t not  too  large) 
and  Eq.  (34)  describes  a field  that  is  relatively  slowly  varying  with  angle. 
The  original  transfer  equation  is  recovered  by  adding  Eqs.  (33)  and  (34). 

Along  with  this  separation,  the  boundary  conditions  are  also  modi- 
fied. In  particular,  a condition  is  placed  on  that  suppresses  back- 
scattering.  That  is,  Lj(0,  p,  0)  = 0 for  p > 0.  Coupled  with  the  condition 
Lj(0,  p,  0)  = 0 for  p < 0 imposed  by  the  fact  that  Lj  is  a diffuse  field  yields 
the  general  result  L^(0,p,  0)  = 0 for  all  p.  This  diffusivity  agreement  also 
requires  L^O.p,  0)  = 0 for  p < 0.  In  addition,  both  Lj  and  L ^ must  remain 
finite  for  t -•  », 

The  method  of  solution  for  L.  is  complicated,  and  is  only  sketched 
out  here.  Standard'  ’ expansions  of  both  L^  and  p are  made.  The  azimuth 
angle  variation  is  handled  by  expansion  in  a series  of  cosine  functions  of 
the  form  cosm0  and  the  zenith  angle  variation  by  expansion  in  a series  of 
associated  Legendre  polynomials  P^(p).  When  these  expansions  are  sub- 
stituted into  the  radiative  transfer  equation  [Eq.  (1)],  a first  order  dif- 
ferential equation  for  the  optical  depth  dependent  expansion  coefficients 

# m 

L,  (t)  can  be  isolated.  This  equation  is  solved  subject  to  the  boundary 
condition  L*  (0)  = 0 for  all  p.  The  solution  is  then  substituted  back  into 
the  expansions  over  i and  m.  The  resulting  form  can  then  be  summed 
explicitly  over  m to  obtain  the  final  solution 

E N u) 

''T/,‘0  £ - 0 <35> 


where  are  the  expansion  coefficients 


\ p(x)P^(x)dx 


for  the  expansion  of  p in  Legendre  polynomials 


(37) 


N 


p(coscp)  = u>xPj(cos<p). 


1=0 

This  contribution  to  the  total  field  L cannot  contribute  directly  to 
the  reflectance  since  L^(0 , p,  0)  = 0 for  all  p.  This  solution  is  required, 
however,  in  the  solution  for  L^.  The  significant  advantage  of  this  approx- 
imation is  that  since  is  assumed  to  be  a slowly  varying  field  with  angle, 
the  solution  of  Eq.  (34)  is  substantially  easier  to  obtain  by  standard  tech- 
niques^^ than  Eq.  (1)  is  for  the  original  field  L.  The  simplest  case, 
that  is,  when  i»  nearly  isotropic,  is  considered  here.  If  is  assumed 
independent  of  the  angle  variables  in  the  integral  source  term  of  Eq.  (34), 
then  the  resulting  transfer  equation  for  is 

dL?(T,p,0)  dL.(T,n,0) 

1 • (38) 


The  solution  of  this  first-order  differential  equation  that  satisfies  the 
boundary  condition  L^O , p,  0)  = 0 for  p < 0 is 


L2(T,p,  0)  = 


(n+n0) 


Lj(t,p,0)  + 


(p+p0Hl-«>0)  E0 


4ir 


N 

X (ZI+IJP^cosVq) 
1=0 


1 

•9 

i 

✓ >n 

*— 

(V  V 

l ) 

— ) 

/ . T TV 

expJ- 

21+1 ) 

V 

\JL1 

ri^a + j.n 

LL  i*  +(- 

1 

1 " 21+1 

i-  i*  i*0-  J 

Finally  then,  the  reflection  coefficient  is  [Eq.  (8)] 

N 


" 4tt[(1  -u»0)p0+p] 


ui 


lPtCOBV  0) 


t-o  ro-«0)i*0 + 0 -T^rW 


(39) 


(40) 


Application  of  this  result  for  p(p,  0)  requires  the  decomposition  of 
the  phase  function  into  a series  of  Legendre  polynomials  according  to 
Eqs.  (36)  and  (37).  This  decomposition  is  considered  to  be  too  complicated 
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for  a model  such  as  we  are  trying  to  formulate  here.  Since  we  have 
assumed  all  along  that  p(<p)  is  a strongly  forward-peaked  function,  we  will 
continue  with  this  assumption  for  the  final  simplification  of  p(p,  0).  Con- 
sider the  analytic  Henyey-Greenstein  phase  function^* ^ 


P('P)  = 


(1-T1  )wQ 
[ 1 - 2r|cos(p  + 


(41) 


whose  expansion  coefficients  are 

U>i  = U)oTli^2jt+1^ 


(42) 


In  the  limit  ri -•  1,  this  function  represents  a strongly  forward-peaked 
function  and  uu^  -♦  u)q(2£+1).  If  this  result  is  used  in  the  denominator  of 
Eq.  (38)  and  Eq.  (39)  used  to  evaluate  the  sum,  then  p(p.,  0)  is  simplified 
to  the  final  result 


p(cos«P0)  p + |xQ 

p(p»0)  p + (l-u>0)p0  ' 


(43) 


The  prediction  of  this  model  for  the  standard  test  case  conditions 
is  shown  in  Fig.  14.  Some  improvement  over  the  simple  single -scattering 
model  is  evident,  but  on  the  whole  the  model  does  not  meet  the  general 
criterion  of  a factor  of  two  accuracy.  Further  discussion  on  this  model 
is  given  in  Section  G. 


F.  Truncation  Approximation 

Hansen^  and  Potter^2^  consider  a straightforward  approximation 
based  on  a simple,  but  physically  justifiable,  redefinition  of  the  problem. 
Consider  a phase  function  p(x)  that  is  strongly  peaked  in  the  forward 
direction  (note,  x = coscp,  <p  = scattering  angle)  and  normalized  by  Eq.  (3) 
to 

1 

C p(x)dx  = 2<»0  . (44) 

-1 


; 
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Separate  off  the  strong  forward  diffraction  peak  by 

p(x)  = u)0A6(l-x)  + p'(x)  (45) 

where  5 is  the  Dirac  delta  function.  Use  of  the  normalization  condition 
Eq.  (44)  gives 

1 

A = j [p(x)  - p‘  (x)  ]dx.  (46) 

0 i 


The  slowly  varying  "remainder"  function  is  normalized  by 


1 

p'(x)dx  = oj0(2-A).  (47) 

-1 

(23) 

The  fundamental  assumption  in  this  approximation  as  put  by  Potter'  ’ 
is  to  "treat  the  radiation  that  has  been  scattered  into  the  forward  direc- 
tion by  the  delta  function  as  not  having  been  scattered  at  all,  " and  to 
still  be  a part  of  the  incident  flux.  Thus,  the  magnitude  and  angular 
distribution  of  scattered  radiation  is  assumed  to  be  controlled  solely  by 
p1  (x).  The  form  of  the  transfer  equation  to  be  solved  is  not  changed.  In 
Eq.  (1),  one  simply  replaces  p by  p'  and  replaces  all  references  to 
by  the  effective  albedo 


u> 


' r,  ai 

= CJD  ’ 1 - — 

0 Oc  2 _ 


(48) 


and,  for  finite  cloud  thickness  problems,  all  references  to  optical  depth 
T by  the  effective  depth 


T' 


(49) 


Since  p'  (x)  is  much  less  rapidly  varying  with  angle  than  is  p(x),  the  solu- 
tion of  the  transformed  problem  is  much  simpler. 
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The  drawback  of  the  method  ia  that  the  truncation  of  the  forward 
diffraction  peak  ia  a subjective  choice.  Moat  applicationa  of  the  method 
have  been  made  in  the  viaible  portion  of  the  apectrum  where  the  forward 
peak  ia  extremely  aharp  and  little  judgment  ia  required  to  determine  the 
beat  truncation.  In  the  infrared,  though  atill  atrongly  peak,  the  truncation 
of  p(x)  ia  ambiguoua.  Becauae  of  thia  ambiguity,  the  method  ia  not  pur- 
sued further. 


G.  Expanaion  in  Ordera  of  Scattering 

The  essence  of  thia  approximation  ia  an  approximate  expansion  of 
the  radiation  field  in  ordera  of  scattering.  Thia  method  has  been  described 
recently  by  Chou*24*  and,  for  the  moat  part,  the  development  here  follows 
that  reference.  Associated  with  each  order  (i.  e.  , once,  twice,  three 
times,  etc.  scattered  photons)  is  an  azimuth  independent  source  function 
Qk(T.p)  and  an  effective  phase  function  p^p,  -pQ>  0Q)  that  is  independent 
of  t.  The  field  is  expanded  as 


Ut,P,0)  = Qk(T'P,Pk*P'0» 


(50) 


k = l 


The  first-order  source  function  and  effective  phase  function  are  obtained 
from  the  single-scattering  solution  of  Eq.  (10)  as 


Q^.p)  =< 


and 


Eor  -T/un  t/u!  Po 

4n  . p + 

p < 0 

Eoe-T/p0  ^0 

4n  P + P0 

b 

p > 0 

Pj  (coscpg)  = p(coa«P0). 

(51) 


(52) 
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The  expansion  of  L(T»Pt0)  by  Eq.  (50)  is  substituted  into  the  transfer 
equation  Eq.  (1)  and  the  following  approximations  made.  First,  Q^(t,ji') 
in  the  integral  source  term  is  assumed  to  be  independent  of  p1,  taken 
outside  the  integral  sign  and  evacuated  at  p'  = p.  Second,  the  effective 
phase  functions  are  defined  by  the  recurrence  relation 

2 it  1 

Pj^H.0#  -KO*0O)  = ^ \ P(n# -Ho,0o)dp'd0'.  (53) 

0 -1 

This  procedure,  along  the  definition 

Qq(t,P)  = E0e”T/^0  (54) 

leads  to  the  following  differential  equation  for  the  source  functions  for 
k 2 0 


dQk+l*T#,A*  1 

•*-*& Qk<T^)=  (55) 


By  working  out  the  solution  for  the  first  few  values  of  k,  one  is  led  by 
induction  to  the  following  solution  that  satisfies  the  boundary  conditions 
Q^(0,  p)  = 0 for  p < 0 and  Q^(T,p)  ->  0 as  t • 


Qk(T,p) 


. E r > _ia_f.-T/w 

° 4tt  p + Pn. 


(56) 


where 


C 

n 


= < 


( 1 ^0 
^ P + PQ 


•) 


n 


p < 0 


(57) 


p > 0. 


♦ 

This  assumption  is  the  major  simplification  of  Chou's  method.  He 
assumes  this  to  be  true  over  the  upward  and  downward  hemispheres  only. 
Here,  the  assumption  of  isotropy  is  made. 
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This  result  gives  the  source  functions  for  use  in  Eq.  (50). 


An  expression  for  p^cos^q)  to  be  used  in  Eq.  (50)  is  obtained  by 
substituting  the  Legendre  polynomial  expansion  of  Eq.  (37)  for  p(cobcPq) 
into  Eq.  (53).  The  result  is 

Pk(c°8V = 4^  (2X+1 ] f4w  22+rf  pj(co‘ V-  < 


The  reflection  coefficient  is  now  obtained  by  substituting  Eqs.  (56)  and  (58) 
into  (50),  evaluating  at  t = 0 for  ji  > 0 and  dividing  by  Eq.  The  result  is 
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Like  the  Romanova  approximation,  application  of  this  model  would  require 
the  expansion  of  p(cos<Pq)  in  a series  of  Legendre  functions.  We  simplify 
the  model  as  in  the  Romanova  approximation  by  use  of  Eq.  (50)  for  in 
the  limit  of  strong  forward  scattering  (ti  "♦  1).  The  result  is,  in  fact,  the 
same  as  that  obtained  in  Eq.  (51)  for  the  Romanova  approximation,  and 
the  consideration  of  accuracy  is  hence  the  same  as  made  in  Section  E. 


5.  SUMMARY 


The  principle  objective  of  this  work  has  been  to  formulate  and  assess 
simple  models  for  computing  the  diffuse  reflectance  of  solar  radiation  from 
plane -parallel  cloud  layers.  By  "simple,"  is  meant  the  ability  to  compute 
the  reflectance  from  no  more  than  the  Mie  scattering  parameters  and  some 
algebraic  formula.  The  most  successful  model  is  the  empirically  modified 
single -scattering  model.  The  nominal  error  of  the  model  is  ± 50  percent. 
Several  other  models  were  considered  (without  empirical  correction)  but 
could  either  not  achieve  suitable  accuracy  in  their  simplest  form  (e.  g. , 
Romanova  and  orders  of  scattering),  were  too  simple  in  their  conception 
(e.  g.  , Turner),  or  were  ambiguous  in  application  (e.  g.  , truncation 
method).  These  results  are  summarized  in  Table  2. 


Table  2 


Summary  of  Reflectance  Models 


Model 

Comments 

Single -Scattering 

Fails  for  u)q  > 0.  5 

Single -Scattering  with  Theoretical 
Modification 

Fails  for  u>q  > 0.  9 

Empirically  Corrected  Single - 
Scattering 

Good  for  all  u)q 

Turner  (two- stream  approximation) 

Fails  for  backscatter 

Order  of  Scattering -Isotropic 

Underestimates  effects  of 
multiple  scattering 

Romanova -Isotropic 

Underestimates  effects  of 
multiple  scattering 

• 

Diffraction  Peak  Truncation 

Manner  of  truncation  is  ambiguous 

* 

I 
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The  Laboratory  Operation*  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 


Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric^chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
l*»ers.  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composite*  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Science*  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetoapheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  opticaL  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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